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(N ■ ABSTRACT 

o , 

Context. Abundance analysis of post-AGB objects as probes of AGB nucleosynthesis. 

Aims. A detailed photospheric abundance study is performed on the carbon-rich post-AGB candidate MACHO 47.2496.8 in the LMC. 
Methods. High-resolution, high signal-to-noise ESO VLT-UVES spectra of MACHO 47.2496.8 are analysed by performing detailed spectrum 
, synthesis modelling using state-of-the-art carbon-rich MARCS atmosphere models. 
O I Results. The spectrum of MACHO 47.2496.8 is not only dominated by bands of carbon bearing molecules, but also by lines of atomic transi- 
' tions of s-process elements. The metallicity of [Fe/H] = - 1 .4 is surprisingly low for a field LMC star. The C/O ratio, however difficult to quantify, 
is greater than 2, and the s-process enrichment is large: the light s-process elements are enhanced by 1.2dex compared to iron ([ls/Fe] = +1.2), 
while for the heavy s-process elements an even stronger enrichment is measured: [hs/Fe] = +2.1. The lead abundance is comparable to the 
[hs/Fe]. With its low intrinsic metallicity and its luminosity at the low end of the carbon star luminosity function, the star represents likely the 
• • ■ final stage of a low initial mass star. 
. ' Conclusions. The LMC RV Tauri star MACHO 47.2496.8 is highly carbon and s-process enriched, and is most probable a genuine post-C(N- 
type) AGB star. This is the first detailed abundance analysis of an extragalactic post-AGB star to date. 

& : 
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1. Introduction 

It is well known that the MACHO experiment towards the 
Large Magellanic Cloud (LMC) yielded invaluable contribu- 
tions to variable star research. The discovery of five distinct 
period-lumi nosity (PL) rela tions of lower mass RGB and AGB 
giants (e.g. lWood et alJl999t) is probably the best known exam- 
ple. Another legacy of the experiment, was the di scovery of 33 
Pop I I Cepheids and RV Tauri stars in the LMC JAlcock et alJ 
Il998h and the definition of a single PL relation of both groups 
with the appearance of the RV Tauri characteristic lightcurves 
(alternating deep and shallow minima) at the high luminos- 
ity end: the RV Tauri stars are a direct extension of Pop II 
Cepheids to longer periods. Absolute luminosities were derived 
and the RV Tauri stars were confirmed to be likely post-AGB 
stars, making these objects members of the small group of post- 
Send offprint requests to: M. Reyniers, maarten@ster.kuleuven.be 
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AGB stars with known luminosities. For Galactic RV Tauri 
stars there is no direct probe t o the luminosity but they were 
identified as post-AGB stars by Jura ( 1986), mainly on the ba- 
sis of the presence of an IR excess due to circumstellar dust 
around many of them. 

The post-AGB phas e is a short phase, so n ot many of these 
sources are known (see IStasiriska et allE00 6. for a catalogue 
of candidate post-AGB objects). Their temperature-gravity do- 
main makes it possible to study a very wide range in chemical 
species by their atomic transitions (if optically bright enough). 
They are therefore ideal objects to constrain the AGB nucle- 
osynthetic and evolutionary models. Interestingly, post-AGB 
stars are chemically much more diverse than theoretically an- 
ticipated: only a very small group of objects shows direct chem- 
ical evidence for AGB nucleosynthesis, being enhanced in car- 
bon and e nhanced in neutron c apture s-process elements (e.g. 
review bv I Van Winckel l2003t and references therein). Since 
the post-AGB evolutionary tracks pass the Cepheid instability 
strip, the pulsating RV Tauri stars could be, in principle, good 
candidates to test further AGB chemical evolutionary models. 
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In recent years it became, however, clear that also for RV 
Tauri stars the chemical picture is complex and chemical ev- 
idence for a post-AGB nature (possible C-enrichement and s- 
process overabundances) is not found in Galactic RV Tauri stars 
(except for maybe V453 Oph, for which a m ild s-process over- 
abundance was found, see iDeroo et alfcooll). In RV Tauri stars, 
depletion abund ance patterns prevail llGiridhar et alJ 120051 
iMaas et aljEool . The basic scenario of the badly understood 
depletion process involves a chemical fractionation due to dust 
formation in the circumstellar environment followed by a de- 
coupling of the gas and the dust with a reaccretion of the 
cleaned gas on the stellar photosphe re, which leav es it de- 
pleted of the refractory elements. Wat ers et al.l ll 19921) showed 
that the most favourable circumstance for this process to occur 
is, if the circumstellar dust is trapped in a disc. The presence 
of a disc in evolved obj ects is likely to be related to binarity 
(e.g. IVan Winckelj l2003h which mea ns that binarity must b e 
widespread in RV Tauri stars (see also lDe Ruvter et all 20061. 

The LMC sample of RV Tauri stars with their known lu- 
minosity is a unique sample to study the nature of these 
stars and of the post-AGB evolution in general. Moreover, 
abundance analyses of post-3rd dredge-up stars in the LMC 
should make it possible to study the yields of the AGB nu- 
cleosynthesis in a more metal deficient environment than the 
Galaxy. With the advent of high-resolution spectrographs on 
8 m class telescopes, it is now feasible to study in detail 
the abundance patterns of these individual RV Tauri candi- 
dates in the LMC. The brightest object of the lAlcock e t al. 
Jl998h sample, MACHO*04:55:43.2-67:51:10, also named 
as MACHO 47.2496.8 , was studied at low resolution by 
IPollard & Llovd Evansl ( 120001) and bv lLlovd Evan7 & Pollard 
( 12004 . The basic parameters of MACHO 47.2496.8 are sum- 
marised in Table ^ Contrary to what is found in Galactic 
RV Tauri stars, this object shows clear indications of a 
strong C-enhancement and of s-process overabundances. The 
low-resolution spectra at different photometric phases of 
MACHO 47.2496. 8 show that it is strongly carbon enriched 
(C/O > 1) with strong C2 bands at deep minima. Spectra taken 
at a resolution of 0.12 nm with the 1.9m SAAO telescope 
around deep minimum (phase 0.88) and with the Anglo- 
Australian Telescope near maximum light (phase 0.61) show 
that the C2 bands disappear almost completely at maximum 
light, with only the actual bandhead of the origin band at 
5 165 A still readily detectable. The spectrum at minimum 
showed strong bandheads of 12 C 12 C at 4737 and 4715 A but the 
4744 A bandhead of 12 C 13 C was very weak or absent, indica- 
tive of a high 12 C/ 13 C ratio. Enhancements of the Ban spectral 
features at 4554 A and 4934 A were detected on these spectra, 
making this object a very interesting star to study the AGB nu- 
cleosynthesis in the LMC in detail. 

In this paper, we will present a detailed abundance analy- 
sis based on high-resolution, high signal-to-noise VLT-UVES 
spectra. The paper is organised as follows: in Sect. |2]we sum- 
marise the available photometry of the object, and discuss its 
variability and extinction. In Sect. [5] the high-resolution obser- 
vations are briefly discussed, while Sect. |4] is devoted to the 
actual analysis and the abundance results. In Sect. [5] we dis- 
cuss the abundance pattern, and address the question whether 



Table 1. Basic parameters of MACHO 47.2496.8. 



MACHO 47.2496.8 


Coordinates 


Q'2000 


4 55 43.23 




02000 


-67 51 10.4 


Mean magnitude 


V 


14.97 


Spectral type 




Supergiant 






pulsating between 






early K and late F 



Source c oordinates: MACHO; magnitud e : lAlcock et all 1 199 8^1: spec- 
tral type: |Pollard & Llovd Evan! <200Cl) . 

the star is intrinsically or extrinsically enriched in s-process el- 
ements. We end with the main conclusions (Sect.[6). 

2. Variability, extinction and luminosity 

2.1. MACHO photometry 

The database of the MACHO project ( lAlcock et all Il992l 
Il995h . accessible via http : //wwwraacho . anu . edu . au/ con- 
tains 260 simultaneous measurements of MACHO 47.2496.8 
in the blue MACHO filter V M and red MACHO filter R M , 
and 465 additional measurements in Vm only. We performed 
a frequency analysis on the 725 Vm magnitudes, using the 
Jurkevich-Stellingwerf phase dispersion minimization (pdm) 
method JStellingweril Il978l) . and deduced a "formal" pe- 
riod (the period between two successive deep minima) of 
P= 1 12.97 d (g-s tatistic of = 0.23) which is slightly longer 
than the period in lAlcock etail Jl998 ). being P = 1 12.45 d. Two 
outliers were removed from the frequency analysis. A part of 
the V M light curve for MACHO 47.24 96. 8 is also shown in 
Fig. 1 of IPollard & Llovd Evansl J2000t) . Our period of 1 12.97 
days and the epoch of minimum light, JD2448673.384, are 
used as the basis for the phases of observations reported here. 

A realistic error estimate on a period determination is al- 
ways difficult to obtain. Here, we use a rather straightforward 
method to give a rough error estimate on the period: the width 
of the pdm minimum is a measure for the spread on the "ac- 
ceptable" periods, i.e. the periods for which the corresponding 
phase diagrams are still acceptable. This width is ~3 10~ 5 c/d 
in frequency space, resulting in an uncertainty of ±0.2 d on 
the period. As an extra check, we performed an additional fre- 
quency analysis on the 260 Vm - Rm colour data points. This 
analysis yielded a slightly shorter period (P=112.86d), but 
still within the quoted error of +0.2 d. 

The conversion of the raw MACHO magnitudes to a stan- 
dard photometric system is not straightforward, since the two 
filters are very broad, and they are only partly overlapping with 
the classical John son or Rron-Cousins filters (see Fig. 1 of 
lAlcock et alll999t for the filter pass bands). The conversion to 
the standard Johnson-Kron-Cousins system is performed using 
the formulae (1), (2), (5) and (6) of lAlcock etail dl999l) . In a 
first step, only the observations that are taken simultaneously in 
both bands, were calibrated to the Johnson-Kron-Cousins sys- 
tem, since Vm - Rm is needed to perform the conversion. Then, 
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Fig. 1. The MACHO V and V-R phase diagram, converted to 
the standard Johnson V and Kron-Cousins R, constructed with 
a period of 112.97 days, and a zero phase at deep minimum 
at JD2448673.384. For the V-R colour diagram, only the si- 
multaneous measurements in both bands are shown. At the 
bottom, the phases of the photometric and spectroscopic data 
are given. The symbols of the spectroscopic data are as fol- 
lows: the full triangles are lo w resolution spectra presented in 
Pollard & Llovd Evans (2000); the open ones are discussed in 
Llovd Evans & Pollard! d2004l) : the phase of the UVES spec- 



trum (0.93) is shown with a hexagon. 

the calibration was extended to all 725 Vm points, using a cal- 
ibration relation that was inferred from the calibration of the 
260 simultaneous points. A phase diagram of the converted V 
magnitude, and the converted V-R colour, is shown in Fig.^ 

From Fig. ^it is clear that the phase diagram displays a 
large dispersion around the mean curve, especially around the 
two minima. This dispersion is caused by cycle-to-cycle vari- 
ations. Note that the colour variations are much more stable 
over the covered cycles. Such cycle-to-cycle behaviour in light 
and colour curves i s often observed in RV Tauri variables (e.g. 
|Pollardetal!l99"a. 



2.2. Near-infared photometry 

Near-infrared photometry was obtained with the Mk III 
Infrared Photometer and chopping secondary at the Cassegrain 
focus of the 1.9 m reflector of the South African Astronomical 



Table 2. Optical Geneva and near-IR photometry of 
MACHO 47.2496.8, together with their JD and phase. The 
Geneva R and I are from the Gunn resp. Cousins photometric 
systems. 



J U 


Phase 


T T R 


\7 

V 


D T 






Geneva 






2453759.6 


0.02 


19.15 16.81 


15.94 


15.18 14.66 


JD 


Phase 


I J 


H 


K 






SAAO 






2450437.4 


0.62 


13.18 


12.67 


12.54 


2450438.3 


0.62 


13.20 


12.56 


12.29 


Mean 


0.62 


13.19 


12.61 


12.41 


<r 




.04 


.04 


.05 






2MASS 






2451111.8 


0.58 


13.12 


12.71 


12.51 






.03 


.04 


.04 






DENIS 






2450377.8 


0.09 


14.02 13.22 




12.53 






.03 .08 




.13 



The Geneva R can be converted to Kg, the Cousins R m agnitude, 
by the conversion forumula from ISchombert etail Jl990h . yielding 
R c = 14.89. 



Observatory by one of us (KRP). This photometry was comple- 
mented by near-IR survey photometry of 2MASS and DENIS. 
All near-IR data, together with their JD and phase, are gathered 
in Table|2l 



2.3. Optical CCD photometry with Euler 

Optical CCD photometry in five different bands was collected 
with the C2 camera at the 1 .2m Euler telescope on La Silla on 
23/24 January, 2006. The filters that were used are the Geneva 
U, B and V, completed with a Gunn-R and a Cousins-I fil- 
ter. Absolute photometry was derived by using exposures of 
the standard stars HD 47645 and SAO 131547 that were taken 
immediately before and after the MACHO 47.2496.8 frames. 
This absolute calibration was cross-checked by the extrac- 
tion of two stars on the LMC-frame with published photom- 
etry which are situated less than 4 ' from MACHO 47.2496.8 : 
CSI-67-045523 (filters U, B and V, Schmi dt-Kaler et al J 19991) 
and MACHO 47.2496. 19 (filters V and R. iKeller et alJl2002l 
Taking into account that there are several transformation for- 
mulae involved in the conversion of the Ge neva and Gunn 
filters to the Johnson-Co usins system ( Schom bert et al] fl990: 
lHarmanec & BozidEoOll) . the magnitudes of the two check- 
stars are very close to the published ones. 

Since there is a nearby star at approx. 2", the extraction of 
MACHO 47.2496.8 was by no means easy. We corrected the 
pixels of our target that are affected by this nearby star by their 
point symmetric counterparts on the unaffected side of the pho- 
tocenter. This induces an extra uncertainty in the derived mag- 
nitudes, but, since the corrected pixels are situated at the bound- 
ary of the psf, this uncertainty turned out to be not larger than 
0.05 magn. The photometry is given in Table [2] The phase of 
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the Euler observations, 0.02, has an estimated error of -0.07, 
due to the quite large time gap between the MACHO photom- 
etry (on which the period is based) and the C2-Euler observa- 
tions. Note that the MACHO photometry is not affected by this 
nearby star, since the latter star has its own MACHO identifi- 
cation: MACH047.2496.21. 



2.4. Spectral Energy Distribution 

The photometry from Table |2] was used to construct a Spectral 
Energy Distribution of MACHO 47.2496.8. The total redden- 
ing E(B-V) was determined by minimising the difference be- 
tween the magnitudes and a reddened model atmosphere. For 
this model atmosphere, we used the parameters that were deter- 
mined in our abundance analysis (Sect. l4~2l . Note that the op- 
tical Geneva photometry and the UVES spectrum from which 
the model parameters are derived, differ less than 0.1 in phase. 
Such a small phase difference is expected to result in only a 
small shift in atmospheric parameters. The result can be found 
in Fig. |21 an E(B-V) = 0.44 was found. Due to the presence 
of the innumerable spectral lines in the model atmosphere, 
the specific passbands of the filters used were taken into ac- 
count (Geneva syst em: http : //obswww . unige . ch/ DENIS: 
iFouque et al 1200(1 SAAO: blasJl973fc . The integration of the 
scaled model combined with an LMC distance of 50kpc, yields 
a luminosity of -5000 L (M bo i - -4.5). 

A specific problem arose for the I magnitude. The Cousins I 
taken with Euler turned out to be suspiciously low. The Cousins 
I was the only magnitude that could not be checked with a star 
on the CCD frame itself (see SectEHl, The DENIS-I magni- 
tude -in gray on Fig. |3- could not be reconciled with the other 
optical photometry as well, possibly attributable to a phase 
shift. We did not take either I magnitude into account in the 
minimisation. 

We made an additional SED with the MACHO photome- 
try, confined to those measurements with a a phase close to 
the UVES phase (difference <0.05), since the model atmo- 
sphere parameters are based on the UVES spectrum. The SED 
made with the MACHO V and R, completed with the DENIS 
I and the near-IR JHK measurements, yielded a reddening of 
E(B-V) = 0.33 and integration of the scaled model yielded a 
total luminosity of -4600 L (M bo i - -4.4). The fit is, how- 
ever, of a lower quality than the fit in Fig. |2] hence we take 
the value of 5000 L as the final value for the luminosity of 
MACHO 47.2496.8. 



2.5. Interstellar extinction 

The reddening derived in Sect. 12.41 is the total redden- 
ing, and can have, besides the interstellar component, a cir- 
cumstellar component as well. The interstellar extinction to- 
wards the LMC is well studied in literature. UBV photom- 
etry and spect ral types of 450 OB st a rs in t he western half 
of the LMC dlsserstedd Il975l 1 19791 1 19821 lArdebere et alJ 
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1972llRousseau et alJl978HCramptorJl979llConti et alll986l 
Fitz patrickl Il988l) give a mean E(B -V) of 0.1 3 mag , using 
the method and calibrations of Fitzpat rick & Gar man 



10' 



10' 



MACH047. 2496.8 

E(B-V) = 0.44 




0.1 



1.0 

Wavelength [/im] 



10.0 



Fig. 2. The spectral energy distribution (SED) of 
MACHO 47.2496.8. Diamonds are the measured magni- 
tudes (converted to fluxes): Geneva U, B, V and Cousins R 
taken with C2+Euler, I from DENIS (in gray), and SAAO 
J, H, K (see also Table 0. The MARCS model is shown in 
gray, while a smoothed version is shown in a full black line. 
The minimisation is made using the unsmoothed model in 
combination with the specific passbands of the photometric 
filters involved. The I magnitude could not be fitted, and is 
excluded in the minimisation. Possibly a phase difference is 
causing this discrepancy. 

There are larger values in OB associations and a region of 
lower values in the North of the LMC. The region surround- 
ing MACHO 47.2496.8 appears fairly uniform on the Digitised 
Sky Survey, and the mean reddeni ng of 16 stars within 30 ar- 
cmin is 0.1 1 mag. However Zari tskv et alJ J2004I) find that the 
distribution of reddening values for cool stars differs from that 
of hot stars and is bimodal, such that stars have E(B-V) near 
zero or near 0.15 mag, depending on whether they lie in front 
of or behind a widespread thin dust layer. We cannot place 
MACHO 47.2496.8 or other carbon stars relative to this dust 
layer and, therefore, have to assume that the reddening is not 
significantly different from the average for carbon stars. 

3. High resolution observations 

High resolution, high signal-to-noise optical spectra of our pro- 
gramme star were taken with the UVES spectrograph mounted 
on the VLT-UT2 (Kueyen) telescope, in the program to ob- 
tain high-quality data of a larger sample of post-AGB objects. 
The data were acquired in visitor mode during ESO period 
#74 by one of us (HVW). These observations fit in the frame- 
work of our ongoing program to study the photospheric chem- 
ical c ompositi on of stars in their last st ages of evolution (e.g. 
iRevniers et alJEoO^ iDeroo et alJl2005h . The resolving power 
of the UVES spectra varies between -60,000 and -65,000. 
Spectra were taken with the chip unbinned. The spectral inter- 
val covered and some other details about the observations are 
given in Table[3] Since there is another star at approx. 2" in the 
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Table 3. Log of the high-resolution VLT-UVES observations. 
Small spectral gaps occur between 5757 A and 5833 A and 
between 8521 A and 8660 A due to the spatial gap between 
the two UVES CCDs. Since continuum spectral intervals are 
hardly found in the spectrum, the signal-to-noise ratios S/N 
given in the rightmost column, are only indicative, and should 
be interpreted more as lower limits than as fixed values. 



date 


UT exp.time wavelength 
start (sec) interval (A) 


S/N 


VLT-UT2 + UVES 


2005-02-08 


02:49 7200 4780-6808 


70 


2005-02-09 


00:30 7200 3758-4983 


40 


2005-02-09 


00:30 7200 6705-10084 


80 



southwest of MACHO 47.2496.8, the slit was placed perpen- 
dicular to this direction. With this instrument setup, the contri- 
bution of the close neighbour is negligible. 

The reduction of our spectra was performed in the ded- 
icated "UVES context" of the MIDAS environment and in- 
cluded bias correction, cosmic hit correction, flat-fielding, 
background correction and sky correction. We used optimal ex- 
traction to convert frames from pixel-pixel to pixel-order space. 
The spectra were normalised by dividing the individual orders 
by a smoothed spline function defined through interactively 
identified continuum points. For a detailed desc ription of the 
reduction procedure, we refer to iRevniersI d2002|) . In Table [3] 
we also list some indicative signal-to-noise values of the final 
data product. Sample spectra can be found in Figs. |4]|5] and [6] 
We note that the most delicate step in the reduction procedure 
is the normalisation, especially in the blue part of the spec- 
trum, since the spectrum is so crowded that the continuum is 
seldom reached. The continuum placement is the most impor- 
tant source of uncertainty on the abundances derived from lines 
in this region. 



4. Abundance analysis 

4.1. Carbon rich model atmospheres 

Before the actual determination of the atmospheric parameters, 
a new grid of models for cool carbon-rich stars was calculated 
with the MARCS code, using opacity sampling in 1 1 000 fre- 
quency points. Atomic, diatomic and polyatomic (C2H2, HCN, 
and C3) opacities were included in addition to the continuous 
ones in the model atmosphere computations. The models are 
hydrostatic, spherically symmetric and computed under the as- 
sumption of LTE; convective energy transport is included us- 
ing a Mixing Length Theory (MLT) formulation. Further de- 
tails on the carbon-st ar model atmospheres will be given in 
I.Tc>rgense n et alJ J 20Q6I): the current grid of models is described 
bv lGustafssonet alJ J2OO3I) . The sub-grid used here was com- 
puted for a rather wide range in T e ff , log g, overall metallicity 
and C/O ratios; the stellar mass was set to 1 M . 



4.2. Atmospheric parameters 

The precise determination of the atmospheric parameters was 
by no means straightforward. The main difficulty was the 
lack of existing comparison spectra. Due to the high C/O ra- 
tio as inferred from the strong C2 bands, the spectrum of 
MACHO 47.2496.8 mimics the spectrum of a cooler carbon 
star in those spectral regions where molecular lines domi- 
nate. The atomic transitions, however, indicate a significantly 
higher temperature than the typical carbon AGB star tempera- 
tures (2500-3500 K). The presence of the molecular lines itself 
points to a temperature lower than the galactic F-type s-process 
enriched post- AGB stars. 

The final atmospheric parameters were derived by an it- 
erative process of fitting specific spectral regions and fine- 
tuning the parameters consecutively. These specific regions 
are carefully chosen in the sense that they contain spec- 
tral features strongly depending on the parameters, like lines 
of different ionisation stages of the same element. Also the 
wings of the Balmer lines Ha and H/? were fitted. The fi- 
nal parameters that were adopted further in the analysis are: 
T eff = 4900 K, log g = 0.0, £ = 3.5 km s _1 and model metallicity 
[M/H] = —1.5. An additional and independent check was made 
by an equivalent width study of selected iron lines. The search 
for clean, unblended lines is not easy in the crowded spectrum 
of MACHO 47.2496.8, but a detailed study of 8 Fe 1 lines and 5 
Fe 11 nicely confirms the atmospheric parameters as derived by 
the spectrum synthesis iteration. 

The procedure of finding the atmospheric parameters of 
MACHO 47.2496 .8 was also followed in previous analys es of 
very similar stars ( Abia et al J20()3lde Lavernv et all200d> . We 
refer to the detailed error analyses in those papers for error es- 
timates on the parameters. Here we only remind that the typical 
errors are as follows: +250 K on T e ff, ±0.5 on \ogg, 1 kms -1 
on the microturbulent velocity £ t and 5 km s _1 on the macrotur- 
bulent velocity. The latter parameter does not influence the line 
strengths, but only the line profiles, and is only used to better 
match synthetic and observed spectra in the spectral syntheses. 

4.3. Spectral regions 

We have analysed in full detail several spectral ranges from 
4000 A to 8000 A. Since the molecular lines of the carbon 
bearing molecules are comparable in strength to those in 
carbon stars, we particularly focussed our spectral synthe- 
ses on those spectral regions that are also used i n the anal- 
ysis of the cooler galactic c arbon stars (see e.g. lAbia et all 
2002i: Ide" Lavernv et alJ l2006). For these regions, line lists are 
carefully compo s ed an d these are discussed in Section 3 of 
de L avernv et alJ (|2006) and will not be repeated here. In gen- 
eral, however, the atomic data of many s-process transitions 
occuring in the whole spectral region covered by our spectra, 
are not known to a good precision, which prevents accurate 
quantitative abundance studies. The most important regions 
are (i) between 4750 A and 4950 A for the s-process elements 
and the mean metallicity, (ii) 6700-6730 A for the Li line, (iii) 
4050-4060 A for Pb, (iv) 4260-4270 A for Tc, (v) the C 2 band 
head around 4737 A for the carbon abundance and isotopic ra- 
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tio 12 C/ 13 C (vi) the CN complex between 8000 and 8100A 
for the N abundance and as an additional tool for the car- 
bon isotopic ratio. The oretical spectra were c omputed wit h the 
TurboSpectrum code (A lvarez & Pledfl998t IPlez et alJ Tl992: 
IPlezI 19921 IPlez et al.ll993l) in spherical geometry. The macro- 
turbulence used ranges between 10 and 20kms~'. 

The results of our abundance analysis are summarised in 
Table 0] and are also graphically presented in Fig. [3] The first 
column in Table 0]gives the solar abundances log e of the el- 
ements that were studied; the second one contains the actual 
ions; in the third one the method of abundance determination 
is given (ss for spectrum synthesis and ew for the equivalent 
width method); the fourth column gives the absolute abun- 
dances derived log e = log(N(el)/N(H))+12; <x ltl , the fifth col- 
umn is the line-to-line scatter of the n lines used (sixth column); 
the last column gives the abundance relative to iron [el/Fe]. 

The solar abundances needed to calculate the [el/Fel values are 

■ I j i j — 1 1 

taken from Grevesse & Sauval ( 1998), except: N (Hibbert et al. 
Il99lh . Mg ( lHolwegerll200lb . La (Lawleret aPbOOll) . Despite 
the fact that there are more recent values for some of the so- 
lar abundances (especially for the solar CNO), we take these 
references to ensure as much as possible consistency with our 
previous analyses, to be able to compare abundance results ad- 
equately. 

Errors on the reported abundances are difficult to assess, 
since very different sources can contribute to the total error, and 
it is often very difficult to perform an adequate error analysis. 
Main sources of error certainly include continuum placement, 
undetected molecular and/or atomic features, an uncertain ef- 
fective temperature and inaccurate \og(gf) values. For error es- 
timates on the reported abundance s, we refer again to the e rror 
analyses in lAbia etal ] J2002I) and lde Lavernv et alJ J2006l) . In 
Fig. [5] typical error bars are drawn, based on the number of 
lines used: if only one line is used, an error of 0.3 dex is ap- 
plied, if more than one line is used, an error of 0.2 dex is as- 
sumed. The latter error estimate is justified if one assumes that 

the error on an abundance based on n lines, goes as -J^of l n - 
In the following subsections, we will briefly discuss our main 
findings for the different elemental groups. 

4.4. CNO 

Apart from the "classical" atmospheric parameters T e « and 
log g, the C/O ratio or, more specifically, the difference between 
the C and O absolute abundances, is by far the most important 
"extra" parameter determining the spectrum of a cool carbon- 
rich star. Unfortunately, it is impossible to determine this differ- 
ence accurately, since the O abundance cannot be determined 
independently; spectral syntheses are very similar for differ- 
ent C and O abundances, as long as the difference e(C)-e(0) 
is kept more or less constant (even within 1 dex of variation 
of both abundances figures). The reason is that extra C and O 
atoms which are added to the photosphere in similar amount, 
immediately combine in the strong CO-bond. Therefore, only 
a lower limit could be derived for the C/O ratio: C/O > 2. The 
carbon isotopic ratio is easier to determine and from the syn- 
thesis of the band head region around 4737 A (Fig|4}, this ratio 



Table 4. The abundance results for MACHO 47.2496.8. The 
explanation of the different columns is given in the text 
(Sect. l4~3l . At the bottom, the usual s-process indices, together 
with metallicity, are given. More information on the calculation 
of these indices is found in the Sect. 14.71 The elements C and 
O are not separately included in the table due to the specific 
difficulties for these elements that are discussed in Sect. 14.41 
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was found to be 12 C/ 13 C = 200 + 25. The error of + 25 indicates 
the formal error due to model parameters (T e g, logg, [M/H]), 
adding the problem of the continuum location and the uncer- 
tainties in the absolute C and O abundances. 

An upper limit for the nitrogen abundance is derived from 
the CN complex at 8000 A, and is found to be N/H < 2.8 10 5 
(log e(A0 < 7.45). The spectrum synthesis of this region pro- 
vides also an additional test for the 12 C/ 13 C ratio. The value of 
200 is compatible with the 12 CN/ 13 CN features in the 8000 A 
range. The derivation of the carbon isotopic ratio is, however, 
much more difficult in this 8000 A region, since the 13 CN fea- 
tures are much weaker than in the 4730 A range, and the uncer- 
tainty of the fit is therefore very large. 
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Fig. 3. The abundances of MACHO 47.2496.8 relative to iron 
[el/Fe]. 



and log e(Fe n) = 6.09. We will refer to the latter value as "the" 
metallicity of MACHO 47.2496.8: [Fe/H] = -1.42. The iron 
abundance derived from the Fe n-lines is chosen, since also at 
the relatively low effective temperature of MACHO 47. 2496. 8, 
the bulk of the iron atoms in the photosphere is singly ionised. 
Note that the Fe proxies nickel and zinc nicely follow the iron 
deficiency. 

The a-elements show surprisingly different abundances: 
Mg seems to be enhanced by 0.3 dex, while Ca is deficient by 
approximately 0.5 dex. The reason is not clear, but the differ- 
ence is likely the initial composition of the object without ex- 
cluding an observational error due to blending: note that we 
derive Mg and Ca abundances using only one line. The lighter 
a-elements, like Mg, are thought to be mainly produced dur- 
ing the hydrostatic burning of massive stars, while the heav- 
ier a-elements come from explosive synthesis during SNe II. 
The abundances of the a-elements in MACHO 47.2496.8 are 
in agreement with the observati onal trends for the LMC that 
were recently derived by PompeiaetajJ ( 120061) . so there is 
no evidence for an intrinsic a-enrichment or depletion. The 
same h olds for Na. Only C a is somewhat lower than ex- 
pected: Pom oeia et al.l (2006) predict [Ca/Fe]~0 for a metal- 
licity of [Fe/H] = -1.4, while we derived [Ca/Fe] = -0.5 for 
MACHO 47.2496.8. 



1 i 1 1 1 i 1 1 1 i 1 1 1 
12 "MACHO 47.2496.8 

. C 2 synthesis (detail) 
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Fig. 4. The VLT-UVES spectrum of MACHO 47.2496. 8 
(points) around the C2 4737 A band head, overplotted with two 
spectrum syntheses with a different 12 C/ 13 C isotopic ratio. 

4.5. Metallicity and a-elements 

The metallicity was already roughly estimated in the iterative 
process of the model atmosphere parameter determination. A 
careful synthesis of four clean Fe 1 lines yields log e(Fe) = 6.0, 
or [Fe/H] = — 1.5. We also calculated the abundance from 8 un- 
blended Fei lines and 5 unblended Fen lines through their 
equivalent width, and the results agree nicely: log e(Fe 1) = 5.98 



4.6. Lithium 

Although the region around the lithium doublet at 6708 A is not 
very well fitted, there is no evidence for an enhanced Li abun- 
dance. The upper limit derived for Li is loge(Li) < 1.5. Note 
that the Li abundance for a strongly s-process enriched object 
is difficult to determine, since the Li doublet is strongly blended 
with aCen transition at 6708.10 A jRevniers et al.l2002l) . 

4.7. S-process elements 

Although trace elements by nature, the absolute dominance of 
transitions of s-process elements in the spectra is overwhelm- 
ing. This is illustrated in Fig|5] in which the synthesis is per- 
formed with and without s-process enhancement. S-process el- 
ements are classically divided into the light ones noted Is (Sr, 
Y and Zr) and the heavy ones, noted hs (Ba, La, Ce, Pr, Nd 
and Sm). Some elements show a wealth of atomic transitions 
in the optical, like Ce and Nd, and the abundances of these el- 
ements can be derived with great accuracy. Others, like Sr and 
Ba, only produce a few strong resonance lines that are not suit- 
able for abundance determination purposes. This fact should be 
taken into account when comparing individual s-process abun- 
dances. We could also derive an abundance for the trans-Ba 
element hafnium (Hf), but it is based on only one line, so it 
should be treated with caution. 

We have also studied two resonance lines of Ban to esti- 
mate the Ba abundance. On top of being strongly saturated, the 
profiles are broader than other atomic lines and both lines show 
a blue-shifted component (14.6 and 15.5 kms -1 , respectively). 
This is probably due to the dynamical structure of the outer 
atmosphere of this pulsating star. Using only the component 
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at the same photospheric velocity as the other atomic lines, we 
derive an upper limit of log e(Ba) < 2.8. We conclude, however, 
that the special line profile makes the Ba abundance very un- 
certain. 

All s-process abundances are based on lines of singly 
ionised species. Additionally, we have also searched for lines 
of doubly ionised Pr and Nd, which are known to produce some 
strong lines in the s pectra of enriched post-AGB stars (e.g. 
iRevniers et ai1l2004 . We found one Ndm line at 5294.099 A 
and two lines of Prm (at 5284.693 A and 5998.930 A) suitable 
for abundance determination. The Ndm line yields the same 
abundance as the one from the Ndn lines reported in Table 0] 
the two Prm lines yield an abundance of log e(Pr) =1.85 and 
2.26 respectively, and are much less consistent with their 
singly ionised counterparts. It is not clear why the discrepancy 
(~0.6dex) is so large. 

In order to characterize the s-process pattern, and to be 
able to adequately compare the observed pattern with other 
s-process enriched stars, three s-process indices are usually 
defined: [ls/Fe], [hs/Fe] and [hs/ls]. Which specific elements 
are taken into account to determine these indices, are, un- 
fortunately, author-dependent, and are mainly just determined 
by the elemental abundances that are obtained in the analy- 
sis. Here, we define the ls-index as the mean of Y and Zr; 
the hs-index as the mean of La, Ce and Nd; and obviously 
[hs/ls] = [hs/Fe] -[ls/Fe]. The elements Pr and Sm are excluded 
from the hs-index to be able to compare with other s-process 
enriched stars. Indeed, in many of the s-process studies, abun- 
dances for these elements are lacking, and hence they are not 
incorporated in the indices. The resulting indices are given at 
the bottom of Table g| 

4.8. Technetium 

The presence of technetium (Tc, Z = 43) is a well known in- 
dicator that a star is currently undergoing dredge-ups on the 
AGB phase, or that it has just left this phase, since the half- 
life of "Tc is ~2x 10 5 yr. The detection of Tc is, however, 
complicated since all available lines are quite weak, and they 
are often blended by much stronger atomic and/or molecular 
lines. One of the lines that is often used in abundance stud- 
ies of AGB stars, is the resonance line at 4262.27 A. The con- 
tinuum placement is the main source of uncertainty in this 
crowded spectral region. Moreover, at the effective tempera- 
ture of MACHO 47.2496.8, the Tc line is intrinsically too weak 
to see a clear difference in the synthesized line for different 
abundances. We have to conclude that Tc cannot be used for 
MACHO 47.2496.8 as a criterion for the intrinsic enrichment. 

4.9. Niobium 

Besides Tc, there is another element whose abundance can be 
used as a criterion to determine whether the star is intrinsi- 
cally or extrinsically enriched: niobium. Nb, which has only 
one stable isotope, is exclusively produced by the decay of 
93 Zr (half-life t i/2 ~ 1.5 x 10 6 yr) and therefore Nb popu- 
lates the stellar surface only million years after the s-process 



elements have been dredged-up. As a consequence, intrinsi- 
cally (post-)AGB stars have a Nb abundance much lower than 
extrinsically enriched stars for which [Nb/Zr] 0. For intrinsi- 
cally enriched stars, the [Nb/Zr] abundance depends only very 
slig htly on the metallicity, and is predicted to be [Nb/Zr] - 1 .0 
(Bist erzo et al"ll2006h . For MACHO 47.2496.8, we checked all 
available Nb i and Nb n lines in the spectrum, but unfortunately 
none of them is suitable to derive a reliable abundance or upper 
limit, since they are all situated in the blue part of the spec- 
trum for which the molecular veiling is most severe. There is 
one Nbn line at 4789.89 A that could give an estimate for the 
upper limit. Again, the continuum position is the main source 
of uncertainty in this spectral region. Also, the blend in the im- 
mediate vicinity of the Nb line at ~4790.5 A is clearly lacking 
in our line list. The derived upper limit of [Nb/Zr] < -0.3 is 
therefore very uncertain, but points rather to an intrinsic than 
an extrinsic nature. Spectra taken in a hotter phase might be 
more suitable to study this Nb line in more detail. 

4.10. Lead 

The end point of the s-process path is situated at lead and 
bismuth. Sinc e the d iscovery of the so-called "lead stars" by 
IVanEcketal1d200ll) . lead has become a very important issue 
in the study of the s-process. Unfortunately, lead has only a few 
lines in the optical region, and they are all situated in crowded 
regions in the blue spectral domain. In MACHO 47.2496.8, 
only the line at 4057.81 A is accessible, and even there the 
blending by molecular lines and s-process elements is severe. 
As a consequence, the main source of uncertainty of the Pb 
synthesis in this region is the continuum placement. Therefore, 
we can only derive an upper limit for the Pb abundance: 
log e(Pb) < 2.8, as illustrated in Fig.|6] 

5. Discussion: intrinsic or extrinsic? 

Our abundance analysis (Table |4] and Fig. [SJl clearly shows 
that MACHO 47.2496.8 is a metal poor, but strongly s-process 
and carbon enriched object in the LMC. Contrary to Galactic 
RV Tauri stars (see Sect. [0, there is no evidence for depletion 
of refractory elements. In a carbo n-rich (C/Q> 1) environm ent, 
ZrC is likely to condense first dLodders & Feglevlfl999l) fol- 
lowed by TiC and SiC. For MACHO 47.2496.8, the abundances 
of these elements follow the trend which is expected from a star 
that underwent an efficient dredge-up while evolving on the 
AGB, so no depletion has taken place. To determine whether 
this enrichment is intrinsic or extrinsic is not straightforward 
but in the following we argue in favour for the intrinsic inter- 
pretation. 

5.1. Abundances 

From Sect. I4.8I it is clear that the spectral feature, Tc, that is 
classically used to determine whether a star is intrinsically or 
extrinsically enriched, cannot be used for MACHO 47.2496.8, 
mainly due to its relatively high T e ff . An alternative indicator, 
the [Nb/Zr] abundance discussed in l4~9l is also difficult to quan- 
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Fig. 5. The VLT-UVES spectrum of MACHO 47.2496.8 (points) overplotted with two spectrum syntheses. Both syntheses are 
made with the same model atmosphere, only the abundances of the s-process elements differ. The synthesis in the dotted line is 
made with solar abundances for the s-process elements (obviously scaled down to the metallicity of MACHO 47.2496.8), while 
the synthesis in the full line is made with the enhanced s-process abundances as tabulated in Table [4] The observed spectrum is 
very well reproduced by the spectrum synthesis, although there are clearly still some lines missing in the line lists. 



tify, but might have some potential for spectra that are taken in 
a hotter phase. 

5.2. Carbon star luminosity function 

Since carbon star counts in the Magellanic Clouds are consid- 
ered to be nearly complete, the luminosity functions for the 
carbon stars (CSLF) in the Clouds can be constructed. In the 
LMC, 7750 carbon stars are detected, and a hi stogram of the 
observational CSLF can be found in e.g. Fig. 6 of Groenewegen 
This histogram does not indicate a clear faint limit to 
the luminosity function of the N-type carbon stars. The re- 
production of the CSLF of the LMC by detailed nucleosyn- 
thetic models is a long standing problem, since these models 
do not succeed in producing the carbon stars at the low lumi- 
nosity end of the CSLF. Recently, it has been argued that the 
low-luminosity end of the CSLF is populated by first giant- 
branch stars that are enriched by mass transfer from a former 
AGB companion. This idea was quantitatively worked out in 



llzzard & Toutl |2004). In this paper, the transition between an 
extrinsically enriched carbon star and a genuine AGB carbon 
star is situated around ~ -4. Consequently, the bolomet- 
ric magnitude M bol =* -4.5 of MACHO 47.2496.8, inferred in 
Sect. |2 does not provide a strong constraint whether the object 
is either intrinsically or extrinsically enriched. 

We compared the flux of our program star to the flux of 
the N stars by integrating over the wavebands B to K. The 
mean figure for MACHO 47.2496.8 has been calculated assum- 
ing that the ratio between maximum and intensity-mean light, 
calculated for V and R using the intensity-scaled light curves 
from the MACHO website, applies to all wavelengths. This is 
a conservative assumption, as the amplitude is probably less in 
the near infrared. The relative flux in BVRJ for the N stars has 
been a ssumed to be th e same as for galactic N stars, using the 
data of Walker Jl979h . Fluxes have been calculated in magni- 
tudes on an arbitrary scale, and show that MACHO 47.2496.8 
is 0.09 mag (at minimum light) and 0.30 mag (at maximum 
light) brighter than the faint limit for N stars in the LMC. 
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Fig. 6. Synthesis of the Pbi line at 4057.81 A. The exact po- 
sition of the continuum is far from clear in the vicinity of the 
line, which makes it only possible to derive an upper limit for 
the Pb abundance. 



It is clear that MACHO 47.2496.8 is located at the low 
luminosity end of the carbon star luminosity function of the 
LMC. Its metallicity is, however, quite low compared to the 
bulk of the LMC stars. A better luminosity function may there- 
fore come from the SMC where it is well known that the CSLF 
is shifted to lower luminosities l lGroen ewegen 1999). The low 
intrinsic metallicity therefore favours even more the intrinsic 
nature of the s-process enrichment. In fact, at the metallicity 
of MACHO 47.2496.8 (~ Z o /30), the minimum luminosity for 
the formation of an AGB carbon star may be slightly lower 

than Mbd 4.0, th e exact value depend ing on the mass-loss 

parametrization (e.g. lStraniero et all2003h . 

5.3. Comparison with galactic CH stars 

CH-stars are metal poor, carbon rich giants with spectral types 
from G to K, with enhanced CH, C2 and s-process elements 
relative to normal giants. CH-stars have typically also large 
radial velocities of ~100kms~ 1 . It is likely that all CH-stars 
are binaries in which the present primary has gained processed 
mate rial from a formerly- AGB c ompanion, now a white dwarf 
(Mc Clure & Woodsworthlll990l) . The first accurate quantita- 
tive abundance analyses were done bv lVantureNl992 J bid) . He 
found that the s-process abundance pattern of the CH-stars is 
characterised by large s-process enhancements and large [hs/ls] 
ratios of >+0.9. Their metal deficiency, carbon enrichment 
and s-process pattern are similar to MACHO 47.2496.8, and 
they are therefore ideal galactic comparison stars. The majority 
of CH-giants, however, show very low carbon isotopic ratios, 
typically 12 C/ 13 C< 10 JVanturelll992al) . which is not seen in 
MACHO 47.2496.8. 



Mor e recent a bunda n ce analy ses o f CH -stars 

J Van Kck etalJ 12001 1 l200/it lAoki et alJ l200ll I2002I) con- 
firmed the older results of IVanturel (11992a) . and focussed on 
the detection of an enhanced lead content, which was predicted 
in different nucleosynthetic AGB models (e.g. iGallino et alJ 
I l998t iGorielv & Mowlavil l2000h . As first discovered by 
I Van Eck et alJ 11200 II) . a subgroup of the CH-stars does show 
this predicted Pb enhancement (typically [Pb/hs] +1), but 
there are also CH-stars that are not compatible with the predic- 
tions, showing [Pb/hs] ratios that are more than one m agnitude 
lower than predicted (see e.g. Fig. 4 in I Van Eck et all2 003). 

The Pb abundance of MACHO 47.2496.8 is very uncer- 
tain, and in Sect. 14.101 we were only able to derive an up- 
per limit of [Pb/hs]~0. Therefore MACHO 47.2496.8 is con- 
sidered as a "lead-low" star. Nevertheless, the low [Pb/hs] ra- 
tio as well as the other s-process indices are still compati- 
ble w ith recent nuc leosynthesis calculations in low-mass AGB 
stars JGalhno & Bisterzol2006l . 

6. Conclusion 

In this paper, we have presented our study of the brightest 
RV Tauri star that was found in the M ACHO experiment: 
MACHO 47.2496.8 JAlcocket al.lll998l) . The star is strongly 
pulsating with a peak-to-peak amplitude in V of ~1.4magn. 
A SED constructed using recent CCD photometry taken with 
the Euler telescope at La Silla, indicates a reddening of E(B- 
V) = 0.44, while the interstellar reddening towards this region 
in the LMC is E(B-V) = 0.11. The intrinsic bolometric mag- 
nitude Mboi calculated from the spectral energy distribution is 
-4.5, placing the star in the low luminosity end of the carbon 
star luminosity function of the LMC. 

Our high-resolution, high signal-to-noise optical UVES 
spectra are unfortunately taken in a cool phase. The main re- 
sults of our detailed abundance analysis performed by a careful 
synthesis of some regions that are already well studied in galac- 
tic carbon stars, are 

- a surprisingly low metallicity: [Fe/H] = —1.4, with no evi- 
dence of depletion 

- a high C/O ratio, certainly larger than 2, and a carbon iso- 
topic ratio 12 C/ 13 C of around 200 

- a large s-process enrichment with [ls/Fe] = +1.2, and 
[hs/Fe] = +2.1 

- a low Pb abundance [Pb/Fe] =s [hs/Fe] 

We used a new grid of state-of-the-art MARCS models for 
cool carbon rich stars, which was especially constructed for 
this analysis. 

Although we cannot totally exclude that the star is extrin- 
sically enriched by a former AGB companion which should 
be seen now as a white dwarf, the enrichment is likely of 
an intrinsic origin. The luminosity, although quite low, is still 
compatible with an advanced phase of a very low-mass star, 
and the pulsations are prototypical of an RV Tauri like object. 
Unfortunately, the star is too hot to use technetium as a criterion 
of an in situ enrichment. Also the alternative element suitable 
as an intrinsic/extrinsic test, niobium, cannot be used for this 
star. New spectra, taken in a hotter phase when the molecular 



M. Reyniers et aL: Abundance analysis of MACHO 47.2496.8 



11 



veiling is absent, could enable a conclusive test concerning the 
intrinsic/extrinsic dilemma. 

With its low luminosity and certainly very low inital metal- 
licity, MACHO 47.2496.8 represents the final evolution of a 
star which must have had an initial mass very close to solar. 
To our knowledge, this is the first detailed chemical analysis of 
a post-AGB star with a known distance and accurate luminos- 
ity estimate. The large carbon enhancement and the very rich 
s-process nucleosynthesis show that also very low-mass objects 
will undergo strong chemical changes during AGB evolution. 
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